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ABSTRACT 



A two-layer spherical a 2 -dynamo model consisting of an inner electrically conducting 
core (magnetic diffusivity \ and radius rj) with a = surrounded by an electrically 
conducting spherical shell (magnetic diffusivity A Q and radius r Q ) with a constant a is 
shown to exhibit oscillatory behavior for values of (3 = \i/\ and ri/r Q relevant to the 
^ ■ solar dynamo. Time-dependent dynamo solutions require rj/r G > 0.55 and (3 < 0(1). 

For the Sun, rj/r Q is about 0.8 and f3 ~ 10~ 3 . The time scale of the oscillations matches 
the 22 year period of the sunspot cycle for Ao = O(10 2 km 2 s _1 ). It is unnecessary 
to hypothesize an aw-dynamo to obtain oscillatory dynamo solutions; an a 2 -dynamo 
suffices provided the spherical shell region of dynamo action lies above a large, less 
magnetically diffusive core, as is the case for the solar dynamo. 

Subject headings: convection — hydrodynamics — instabilities — magnetic fields - 
Sun: magnetic fields — stars: magnetic fields 



1. Introduction 



a 2 - and aw-dynamo models have been studied for decades (Braginsky 1964; Steenbeck and 
Krause 1966; Roberts 1972; Moffatt 1978; Gubbins and Roberts 1987; Baryshnikova and Shukurov 
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1987). The a 2 -dynamo is usually stationary (e.g., Roberts 1972; Gubbins and Roberts 1987; Holler- 
bach 1996), although oscillatory a 2 -dynamos have been found to occur in the special circumstance 
wherein a changes rapidly in boundary layers (Radler and Brauer 1987; Baryshnikova and Shukurov 
1987). In such cases, the period of the a 2 -dynamo depends strongly on the location of the a- 
boundary layer and is typically an order of magnitude or more smaller than the magnetic diffusion 
time across the dynamo generation region. In general, oscillatory dynamo behavior has been pro- 
duced by combination of the a- and w-effects. Kinematic models of the solar dynamo, which is 
inherently oscillatory, have been of the aw-type (e.g., Roald and Thomas 1997). 

In this paper, we report that spherical oscillatory a 2 -dynamos can be simply induced by the 
magnetic coupling between an electrically conducting outer fluid shell and a conducting inner 
spherical core even when a in the outer shell is a constant. The period of oscillation is of the same 
order of magnitude as the magnetic diffusion time across the outer shell and depends largely on 
the electrical conductivity of the inner core. Oscillatory behavior occurs when the outer region of 
dynamo action surrounds a large, less magnetically diffusive core. The radiative interior of the Sun 
is a large region with a smaller magnetic diffusivity than the overlying convection zone wherein 
dynamo action occurs. The oscillatory character of the Sun's magnetic field, as expressed in the 
22 year periodicity of the sunspot cycle, could then be related to the electromagnetic coupling of 
the region of magnetic field generation in the convection zone with the radiative core of the Sun. 
The importance of an inner electrically conducting core to the problem of magnetic field generation 
in an overlying spherical shell has been emphasized in the a 2 -type models of the geodynamo by 
Hollerbach and Jones (1993, 1995), Hollerbach (1996), and Gubbins (1999). Nevertheless, the effects 
of an inner core, with magnetic diffusivity different from that of the overlying fluid convecting shell 
in which dynamo action takes place, have not been fully elucidated. 

The problem of inner core-fluid shell coupling is made difficult by complicated electromagnetic 
matching conditions at the interface between the regions. For this reason, and to facilitate under- 
standing of the physical effects, we consider the simplest type of a 2 -dynamo model consisting of a 
spherical shell with a = constant surrounding a core with a = 0. We derive the appropriate match- 
ing conditions for a core and shell of arbitrary magnetic diffusivity. These matching conditions do 
not appear to have been considered in previous studies of spherical a 2 -dynamos and they result in 
oscillatory dynamo solutions. 

2. Model, Equations and Boundary Conditions 

The model consists of a turbulent fluid spherical shell of inner radius ri and outer radius r a with 
constant (turbulent) magnetic diffusivity A G . A magnetic field is generated in the shell by the a- 
effect (Steenbeck and Krause 1966; Roberts 1972). For r > r Q , we assume there is a non-conductor; 
for r < Ti we assume that there is a conductor with magnetic diffusivity \. The kinematics of 
the a 2 -dynamo in the spherical shell is governed by the non-dimensional linear equations for the 
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magnetic field B G 

^ = K(l-r,)VxaB + V 2 B (1) 

V • B D = (2) 
In the inner sphere the magnetic field Bj is governed by 

& -§ = /?V 2 B, (3) 

V • B, = (4) 

Equations (l)-(4) are scaled by the thickness of the shell (r G — n) and by the magnetic diffusion 
timescale (r — rj) 2 /A . The scaling of the linear system of equations for the magnetic field is 
arbitrary. The non-dimensional parameters in the above equations, (3, rj, and the magnetic Reynolds 
number R are defined as 

fit Xi n D ^ /r\ 

P = y, V=-, R=~r- (5) 

Since the main purpose of this paper is to understand the effect of an electrically conducting inner 
core, we adopt the simplest possible model and take a constant in the spherical shell r% < r < r a ; a 
is zero outside the shell. With this assumption, spherical harmonics are decoupled and the problem 
is reduced to a one-dimensional problem with complicated boundary conditions. 

At the interface between the shell and the perfectly insulating exterior, i.e., at r = r Q , the 
magnetic field must be continuous 

B = B (e) at r = r Q (6) 

where B^ e ) = — V<p is the magnetic field in the insulating exterior r > r Q , and V 2 ^ = 0. At the 
interface between the shell and the conducting inner sphere, i.e., at r = rj, both the magnetic field 
B and tangential components of the electric field E must be continuous 

B = Bj, f x E = r x Ej at r = n (7) 

where r is the unit radial vector, E Q is the electric field in the outer shell, and Ej is the electric 
field in the inner core. 

Conditions (2) and (4) allow us to express the magnetic fields as a sum of poloidal and toroidal 
vectors 

B = V x V x r/i + V x rg Q (8) 
Bi = V x V x r/ij + V xrgi (9) 

where r is the position vector. Use of Equation (8) in boundary condition (6) and expansion of h 
and g Q in terms of spherical harmonics give 

dh Q (l + l)h 

9o = 0, -5- + = at r = r D (10) 

or r 
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where I is the degree of the spherical harmonic Y'( n . 

Extra care must be taken for the magnetic boundary conditions at the interface r = Tj. There 
are four different cases that we have studied: 

(I) . The limit j3 — > oo for both stationary and oscillatory dynamos. In this case, the boundary 
condition for the magnetic field is simply 

go = 0, ^£-^£ = atr = r l (11) 
or r 

(II) . The limit (3 — > for a stationary dynamo. In this case, boundary conditions (7) require 

h o = 0, ^(1-^)^-^1=0 at r = n (12) 

(III) . The limit (3 — > for an oscillatory dynamo. In this case, boundary conditions (7) require 

^-^ = 0, i?(l-,)(/ + l)/ lo -%l = atr = r, (13) 

(IV) . The general case for (3 not tending toward or oo, for both stationary and oscillatory dynamos. 
In this case, boundary conditions (7) require 

9 o = 9 „ K = h, f= = ^, *(!-,) ^-^ +/3 %»>=0 at r = n (14) 
or Or or or Or 

The last case is evidently the most complicated one. The solutions presented below show that 
for (3 > 10, case I provides a good approximation to case IV, while for (3 < 0.1, cases II and III 
provide a good approximation to case IV. In cases II through IV, g and h are coupled by boundary 
conditions (7). The solutions are invariant to a change in the sign of R. 



3. Solution Method 

In all cases, solutions are expanded in terms of spherical harmonics, implicit in the forms of 
the boundary and interface conditions given above. The spherical harmonics are decoupled and 
only the lowest one (£ = 1) is used in the analysis. The i = 2 mode, not discussed here, behaves 
similarly to the £ = 1 mode. 

The time dependence of the solutions is written as exp (a r + iujt) and onset of dynamo action 
(cr r = 0) is sought. As discussed below, dynamos are either stationary (to = 0) or oscillatory (w^0). 
The frequency u> is dimensionless with respect to the timescale (r — ri) 2 /A . 

In case I, a perfectly insulating core, and in cases II and III, a perfectly conducting core, it 
is only necessary to solve for g Q and h Q subject to the above boundary conditions at and r Q . 
Exact analytic solutions for g D and h D in these cases can be found in terms of the spherical Bessel 
functions of the first and second kind. The solutions reduce to finding the eigenvalues of a 4 x 4 



- 5 - 



matrix. The eigenvalues give critical values of the magnetic Reynolds number R as a function of 
T) = Ti/r a , for which steady or oscillatory dynamos are possible (i.e., a r = 0). 

In case IV, a core of arbitrary /?, solutions must be obtained in both the shell and the core 
subject to the above matching conditions, i.e., g , h Q , gi, and hi must be determined. In prin- 
ciple, analytic solutions are possible, but it is computationally more efficient to seek numerical 
solutions. We do this by employing a spectral- Tau method which expands solutions in terms of 
Chebyshev polynomials. The numerical solutions for arbitrary (3, and the analytic solutions deter- 
mined independently in cases I, II, and III, provide a mutual validation of the separate methods. 
For appropriate values of (3, the solutions of the separate methods agree essentially exactly. 



The principal results of this study are summarized in Figure 1 which gives the critical magnetic 
Reynolds number R cr for the onset of dynamo action in the dipole (/ = 1) mode as a function of 
r), the ratio of the inner radius of the shell to its outer radius. The critical value of R increases 
with increasing rj. For an insulating core ((3 >> 1, dashed curve) dynamo solutions are always 
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Fig. 1. — Critical magnetic Reynolds number R cr and dimensionless frequency lo of dipolar dynamo 
solutions at the onset of dynamo action vs. r]. 
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steady (ui = 0). For a perfectly conducting core ((3 << 1, solid curve) there are two branches of 
dynamo solutions depending on r\\ for r\ < 0.55 the dynamo is steady, but for r] > 0.55 the dynamo 
is oscillatory. There is a jump in R cr at the transition from steady dynamos to oscillatory dynamos 
near r\ = 0.55. The dimensionless frequency ui of the oscillatory dynamos, also shown in Figure 1 
as a function of rj, varies between about 2.5 and 3 for all values of rj considered. 

The values of R cr for arbitrary (3 lie in the narrow space between the solid and dashed curves 
of Figure 1. Importantly, it is found that j3 need not in fact be very small compared with unity 
for oscillatory dynamos to exist. For example, when 77 = 0.8, the value appropriate to the solar 
dynamo, oscillatory dynamo solutions are found for (3 less than 2 to 3. 

5. Discussion 

The problem solved above is a classically simple one of the type considered by Steenbeck and 
Krause (1966) and Roberts (1972) decades ago. Yet the effects of an inner electrically conducting 
core on the a 2 -dynamo are subtle and not heretofore appreciated. They enter through the compli- 
cated electromagnetic matching conditions at the interface between the core and the surrounding 
shell in which dynamo action occurs. The main effect of the core is to introduce time dependence 
into the dynamo solutions for cores whose radii are greater than about 0.55 of the outer radius of 
the shell. An additional requirement for time dependence is that the core be a reasonably good 
electrical conductor; in terms of the magnetic diffusivity ratio (3 = Aj/A , (3 < 0(1) suffices for 
oscillatory dynamo behavior. 

The importance of all this to the solar dynamo is that the parameters of the solar dynamo 
satisfy the requirements of oscillatory a 2 -dynamo solutions. For the solar dynamo 77 is about 
0.8 and (3 is about 10 -3 (Moffatt 1978). In addition, if uj (from Figure 1) is made dimensional 
using the time scale (r D — r-j) 2 /A D with r Q — = 1.4 x 10 5 km and A Q = O (lO 2 km 2 s _1 ) (eddy 
magnetic diffusivity), then the period of the oscillatory dynamo solution is comparable to the 
22 year period of the sunspot cycle. Thus, a 2 -dynamo action alone could be responsible for the 
observed time dependence of the large scale solar magnetic field. It is not our intent to suggest that 
the w-effect is not significant in dynamo action in general, or in the solar dynamo in particular, 
because it represents a physically important process. Our purpose is only to clarify some physics 
and demonstrate the potential importance of a hitherto overlooked effect, that of the oscillatory 
a 2 -dynamo. A detailed analysis of the cases a = a(r) and a cx cos 9(9 = polar angle) is in progress. 

G. S. acknowledges support from NASA's Planetary Atmospheres Program. K. Z. is supported 
by PPARC and NATO grants. 
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